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Abstract. We present a computer-aided, yet fully rigorous, proof of Ira Gessel's tantalizingly 
simply-stated conjecture that the number of ways of walking 2n steps in the region x + y > 
0; V > of the square-lattice with unit steps in the east, west, north, and south directions, that 
start and end at the origin, equals 16™ ^5/3" ' 



1. Introduction 

There is a certain family of lattice walks, let's call them the Gessel walks, whose counting function is 
puzzling the combinatorialists already for several years. Gessel walks (that are trivially equivalent 
to the walks described in the abstract) are walks in the lattice H? that stay entirely in the first 
quadrant (viz. they are walks in N 2 ) and that only consist of unit steps chosen from G := {<—, — 

/*}. If f(n; i,j) denotes the number of Gessel walks with exactly n steps starting at the origin 
(0,0) and ending at the point (i,j), then the counting function is the multivariate power series 

00 00 00 

F(t;x,y) ^EEE^^'W"- 

n=0 i=0 j=0 

We would call this power series holonomic (with respect to t) if it satisfies an ordinary linear 
differential equation (with respect to t) with polynomial coefficients in t, x, y. This may or may 
not be the case. 

For example, the Kreweras walks are defined just like the Gessel walks, but with the unit steps 
chosen from {<— , |, /*} instead of from G. It is a classical result [7] that their counting function 
is holonomic. In constrast, Bousquet-Melou and Petkovsek showed that the counting function for 
certain Knight walks is not holonomic [3] . Mishna [8] provides a systematic study of all the possible 
walks in the quarter plane with steps chosen from any step set S C {<— , "\, f, — \,, J., ,/} with 
\S\ = 3. She shows that the counting functions for the step sets {/*, \, \} and {/*, \, f } (and 
some others that are equivalent to those by symmetry) are not holonomic while all others are 
holonomic. 

For the number of walks returning to the origin, there is sometimes a nice closed form represen- 
tation, even if there is no such representation for the number of walks to an arbitrary point 
For instance, if k(n; i, j) denotes the number of Kreweras walks of length n from (0, 0) to (i, j) [9, 
A006335], then [7] 

4™ /3rA 
fc(3n;0,0) = 7 — (n > 0) 



(n + l)(2n + l) 
and k(n; 0, 0) = if n is not a multiple of 3. 

Gessel [5] observed that a similar representation seems to exist for the number /(n; 0, 0) of Gessel 
walks returning to the origin [9, A135404]. He conjectured the following closed form representation. 

Theorem. Let f(n;i,j) denote the number of Gessel walks going in n steps from (0,0) to 
Then /(n;0,0) = 0ifnis odd and 

/(2 "<°' 0) = 16 (5/3) n (2)„ ( ^ 0) ' 
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where (a) n := a(a + 1) ■ • • (a + n — 1) denotes the Pochhammer symbol. 

The purpose of the present article is to describe, to human beings, the proof of this theorem. The 
proof is accomplished by computing a homogeneous linear recurrence in n for f(n; 0, 0). Then the 
statement follows directly by verifying that the right hand side satisfies the same recurrence and 
that the initial values match. Our recurrence has order 32, polynomial coefficients of degree 172, 
and involves integers with up to 385 decimal digits. As this is somewhat too much to be printed 
here (it would cover about 250 pages), we provide it electronically at 

http://www.math.rutgers.edu/~zeilberg/tokhniot/Guesscl2 

which has the recurrence as a Maple procedure whose output being proves Gessel's conjecture 
(once the first 32 initial values are checked, but this has already been done by Gessel himself, 
when he formuated his conjecture) . 

Our result implies that F(t; 0, 0) is holonomic with respect to t, but it has no direct implications 
concerning the holonomy of F(t;x,y) for other x,y of interest, e.g., x = y = 1, and all the more 
so for general x, y. There is, however, strong evidence that even the general counting function 
F(t;x,y) with "symbolic" x and y is holonomic in t, see [2]. 

2. The Quasi-Holonomic Ansatz 

2.1. Annihilating Operators. Let S n ,Si,Sj be the shift operators, acting on f(n;i,j) in the 
natural way, e.g., S n f{n; i,j) = f(n+l;i,j). An annihilating operator of f(n; i,j) is an operator R 
with 

R(n,i,j,S n ,Si,Sj)f(n;i,j) = 0. 
Those operators belong to a noncommutative polynomial algebra Q(n, i, j)[S n , Si, Sj], and to- 
gether they form a left ideal in that algebra, called the annihilator of f(n; Note that for an 
annihilating operator we can always assume polynomial coefficients instead of rational coefficients, 
i.e., that the operator belongs to (Q[n, i,j][S„, Si, Sj}, by clearing denominators. 
Our goal is to find an annihilating operator for Gessel's f(n; i,j) that implies the conjecture. For 
example, it would be sufficient to know an annihilating operator R(n,i,j,S n ) free of the shifts Si 
and Sj, because then R(n, 0, 0, S n ) would be an annihilating operator for f(n; 0, 0). 
For this reasoning to apply, we could actually be less restrictive and allow also shifts in i and j to 
occur in R, as long as they disappear when i and j are set to zero. Our goal, therefore, is to find 
operators P, Q\, Qi such that 

R{n,i,j,S n ,S t ,Sj) = P(n,S n ) + iQx{n,i, j, S n , Si, Sj) + jQ 2 {n,i, j, S n , S u Sj) 
annihilates f(n;i,j). This is the quasi-holonomic ansatz [6]. 

2.2. Discovering annihilating operators. We search for annihilating operators by making, for 
some fixed d, an ansatz 

R= E c ei _ e6 n*H e >f>S^S?Sf 

0<ei ,...,es<d 

with undetermined coefficients c ei; ... ie6 . Applying this "operator template" to f(n;i,j) gives 

^2 c eu ..., e6 n ei i e2 f 3 f(n + e 4 ;i + e 5 ,j + e 6 ), 

0<ei ,....C6<d 

which, when equated to zero for any specific choice of n,i,j yields a linear constraint for the 
undetermined coefficients. (Note that f(n; i,j) can be computed efficiently for any given n, i, j e 
7L.) 

By taking several different n,i,j we obtain a linear system of equations. If that system has no 
solution, then there is definitely no annihilating operator matching the template. If there are 
solutions, then these are candidates for annihilating operators. 

We can clearly restrict the search to quasi-holonomic operators by leaving out unwanted terms in 
the ansatz for R. 
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2.3. Verifying conjectured annihilating operators. An algorithm was given in [6] for deciding 
whether some given operator R G Q(n, i, j) [S n , Si, Sj] annihilates f(n; i, j) or not. We repeat this 
algorithm for the sake of self-containedness. 

First note that the step set {<—, — >, ,/} gives readily rise to the recurrence 

f(n + = f(n;i + + f(n;i- + f{n;i + 1, j + 1) + f(n;i - l,j - 1), 

and therefore the "trivial operator" 

T . — SfiSi S ' j S ^ Sj Sj S j 1 

certainly annihilates f(n; Instead of checking that R annihilates f(n; we will check that 
TR annihilates f(n;i,j). By the following lemma, this is sufficient. 

Lemma. Suppose that an operator R is such that (T R) f (n; i , j) = 0. Then it can be checked 
whether Rf(n;i,j) = 0. 

Proof. (T R) f (n; i , j) = implies that T annihilates Rf(n;i,j), i.e., Rf(n;i,j) also satisfies the 
above recurrence. Therefore, in order to show that Rf(n;i,j) = entirely, it suffices to show that 
Rf(n; i,j) = for n = and all i and j. If r n bounds the degree of S n in R, then it suffices to 
verify Rf(n;i,j) — for n — and < i,j < r n , because we clearly have f(n;i,j) — for i > n 
or j > n. This leaves us with checking finitely many values, which can be done. □ 

By the lemma, in order to check Rf(n;i,j) = 0, it suffices to be able to check (T R) f (n; i , j) = 0. 
For checking the latter, compute operators U, V with TR = UT + V by division with remainder. 
Then 

(TR)f(n;i,j) = Vf(n;i,j) = 0. 

If V is the zero operator, then we are done, otherwise we proceed recursively to show that 
Vf(n;i,j) = (compute U',V with TV = U'T + V, observe that (TV)f(n;i,j) = iff 
V'f(n;i,j) — 0, and so on.) As T has constant coefficients and, for any d > 0, the commutation 
rules in Q(n,i, j)[S n , Si, Sj] are such that S t n d = n d S l , S-jn d = n d Sj and S n n d = n d S n + 0(n d_1 ) 
(and similarly for i and j in place of n), it follows that the degree of V with respect to n, i,j will 
be strictly smaller than the degree of R with respect to these variables. Therefore, the recursion 
must eventually come to an end. 

2.4. Nice idea, but. . . At this point, we know that all we need for proving the conjecture is a 
quasi-holonomic annihilating operator for f(n;i,j). We know how to search for such operators, 
and once empirically discovered, we know how to verify them. 

Unfortunately, it turned out that if a quasi-holonomic annihilating operator for f(n;i,j) exists at 
all, then it must be quite large. It was shown [6] that there is no such operator of order up to 8 
in either direction with polynomial coefficients of total degree at most 6. Increasing the bounds 
on order and degree further might, of course, help, but this is beyond our current computing 
capabilities. (For the above assertion, a dense linear system with several thousand variables and 
equations had to be solved exactly.) 

3. A Takayama-Style Approach 

By making an ansatz, we could not find a quasi-holonomic annihilating operator, but we could 
find (and verify) plenty of other operators, R\,R 2 ,Rz, . . . that were not of the quasi-holonomic 
type. Once it has been verified that these Ri are indeed annihilating operators, we may of course 
freely choose any operators Pi, P2, P3, ■ ■ ■ , and the combined operator 

Pii?i + P 2 R 2 + P3R3 + ■■■ 

will again be an annihilating operator. In other words, all annihilating operators form a left 
ideal in the corresponding algebra. Our next step is to find a quasiholonomic combination of the 
operators Ri, R2, R3, ■ ■ ■ that were found (and verified) by the method of the previous section. 
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3.1. Takayama's Algorithm. Assume we want to find a recurrence for the sum 

k 

In his "holonomic systems approach" [12] the third-named author proposes to search for an anni- 
hilating operator R of f(n, k) of the form 

R(n, S k ,S n ) = P(n, S n ) + (S k - l)Q(n, S k , S n ). 

Summing over k shows (in case of natural boundaries which we will assume in the following) 
that P(n, S n ) annihilates the sum. Starting with the annihilator of the summand f(n, k) in 
A = Q(fc, n)[Sk, S n ], i.e. Ann^/ C A, one computes an R(n, S k , S n ) £ Ann^ / free of k (e.g., by 
elimination via Grobncr bases). Any such R can be brought to the desired form as above. 

Almkvist [1] observed that in the above setting the constraint for Q can be released: The whole 
proof would go through in the same way if additionally Q depends on k. This fact is exploited 
in Takayama's algorithm [10, 11] which originally was formulated only in the context of the Weyl 
algebra. Chyzak and Salvy [4] extended the algorithm to more general Ore algebras (which include 
also the shift case that we are dealing with) and proposed some optimizations. The idea in short 
is the following: While in the algorithm above, first k was eliminated and then the part (Sk — 1)<3 
was removed (which corresponds to divide out the right ideal (Sk — 1) A), the order is now reversed. 
In Takayama's algorithm we first reduce modulo (Sk — 1)^4 and then perform the elimination of k. 
The algorithm usually leads to shorter recurrences since we allow more freedom for Q. Second, 
the elimination is in general much faster since we got rid of Q from the very beginning. Note that 
Q is not computed at all so we have to assure natural boundaries a priori. 

There is one technical complication in this approach. The fact that we are computing in a non- 
commutative algebra restricts us in the computations after having divided out the right ideal 
(Sk — 1)A In particular, we are no longer allowed to multiply by k from the left. We can easily 
convince ourselves that otherwise we would get wrong results: Assume we have written an oper- 
ator already in the form P + (Sk — 1)Q. Multiplying it by k and then reducing it by (Sk — 1)A 
leads to kP — Q since we have to rewrite k(S k — 1) as (Sk — l)(k — 1) — 1. Because k does not 
commute with S k — 1 we get the additional term — Q in the result which we lose if we first remove 
(Sk — 1)<3 and then multiply by k. 

In order to find a fc-free operator one needs an elimination procedure that avoids multiplying 
by k. Let now Ri, . . . , R m £ A be the operators which generate Ann^ /, and let R[, . . . , R' m £ 
Q(k,n)[S n ] be the corresponding reductions modulo (Sk — 1)A For 1 < i < m we can write 

R'i(k,n, S n ) = Rifi(n, S n ) + Ri,i(n, S n )k + Rip(n, S n )k 2 + . . . 

where Ri j £ A' := Q(n)[5„]. Elimination of k now amounts to finding a linear combination 

/ Ri,o \ 
Ri,i 

Rl,2 



Pi(n,S n ) 



V 



+ ■ ■ ■ + P m (n, S n ) 



) 



( Rm.O \ 


/ P(n, S„) \ 


Rm.l 





Rm,2 





V : / 


V / 



for some Pi,... ,P m £ A. The vector on the right hand side corresponds to the desired k-ivee 
operator. In general, this will not work yet since we can not expect to succeed in the elimination 
without multiplying by k at all. Hence we also have to include multiples of the Ri by powers of k. 
More algebraically speaking, the computations take place in an A-module M that is generated by 
the above vectors plus all elements mod (Sk — 1)A 1 < i < J = 1,2,.... The elimination 

is achieved by computing a Grobner basis of this module. Note that P £ M if and only if there 
exists a Q £ A such that P + (Sk — 1)Q £ Ann^ /. For practical purposes we have to truncate 
the module M by considering only elements up to a certain length d, i.e., which have zeros in all 
positions greater than d. The most natural choice for d is the highest power k d that appears in 
Ri, . . . , R m . But we are not guaranteed that for any P,Q £ A with P + (Sk — 1)Q £ Ann^ / the 
operator P is an element of the truncated module. In the unlucky case that no k-iiee operator 
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was found, the bound d has to be increased. The algorithm works similar in the case of multiple 
sums where we want to eliminate several variables k\, . . . , k r . 

3.2. Proof of Gessel's conjecture. Now back to Gessel's conjecture: Recall that we were look- 
ing for a quasi-holonomic operator 

R(n,i,j,S n ,Si,Sj) = P(n,S n ) + iQi{n,i,j, S n , Si, Sj) + jQ 2 {n,i,j, S n , Si, Sj) 

where we are mainly interested in P, because Qi and Q2 anyway vanish when we set i and j to 
0. This task is very similar to the setting in the previous section and with slight modifications we 
can apply Takayama's algorithm to solve it. The only difference is that now i and j play the role 
of Sk — 1, and instead of k, we want to eliminate the operators Si and Sj. Consequently we have 
to consider the Q(n)[iS n ]-module which is generated by {S^ 1 Sj 2 \e\, ey. = 0, 1, . . .}. 

For our concrete application, we started with a set of 16 annihilating operators for f(n; i, j). These 
operators were found by the ansatz described in section 2.2 and verified as proposed in section 2.3. 
The maximal degree w.r.t. i as well as the maximal degree w.r.t. j is 4. Some of the operators 
had degree less than 4 in i or j, hence we had to add their corresponding multiples to the set of 
annihilating operators (which after this step consisted of 24 elements). Next we performed the 
substitution i = and j = 0. Finally, the elimination of Si and Sj using Grobner bases took 
about 30 hours and resulted in an operator P(n, S n ) of order 32 and polynomial coefficients of 
degree 172 in n. (This is the operator posted on our website.) 

As already pointed out, Takayama's algorithm does not deliver Q\ and Qi. However, in principle 
the full certificate R(n,i, j, S n , Si, Sj) can be computed by doing some book-keeping during the 
run of the algorithm. But in the case of Gessel's conjecture this extra cost would make our 
computations not feasible with current computers (we would have to wait for a few more Moore- 
doublings). We want to emphasize that, nevertheless, the proof is completely rigorous. 

To cite a simple (commutative) analogy, Euclid devised, more than 2300 years ago, an algorithm 
to find the greatest common divisor of two integers. Later mathematicians extended it to the 
generalized Euclidean algorithm that inputs integers m and n and outputs not only d, the greatest 
common divisor of m and n, but also two other integers a and b such that am + bn = d. Just 
because our computers are not fast or big enough to actually find these two other integers a and 
b does not detract from the correctness of the output d of the original Euclidean algorithm, and 
their existence is implied by it. 



4. Conclusion 

Computer-aided proofs have come a long way since the hostile reception of the Appcl-Hakcn 
landmark proof of the Four-Color Conjecture. Even as recently as 1998, Hales' breakthrough 
computer-aided proof of Kepler's conjecture was met with skepticism, but it did get published, 
with some reservations, in the prestigious journal Annals of Mathematics. Both the Appel-Haken 
and Hales theorems are examples of extremely simply-stated statements, whose proofs defied, so 
far, human attepts. While Ira Gessel's conjecture has neither the longevity nor the notoriety of 
the above theorems, it does belong to that genre, and we believe that it is very possible that 
a short human proof does not exist. Unfortunately, formally proving this last meta-conjecture 
would be probably much more difficult than proving Gessel's conjecture, since proving realistic 
lower bounds is a notoriously difficult task. So we have to resort to empirical sociological testing, 
using the ingrained greed of human mathematicians. To that end, the third-named author (DZ) 
offers a prize of one hundred (100) US-dollars for a short, self-contained, human-generated (and 
computer-free) proof of Gessel's conjecture, not to exceed five standard pages typed in standard 
font. The longer that prize would remain unclaimed, the more (empirical) evidence we would have 
that a proof of Gessel's conjecture is indeed beyond the scope of humankind. 



6 



MANUEL KAUERS AND CHRISTOPH KOUTSCHAN AND DORON ZEILBERGER 



References 

[1] Gcrt Almkvist and Doron Zeilberger. The method of differentiating under the integral sign. Journal of Symbolic 

Computation, 11(6):571-591, 1990. 
[2] Alin Bostan and Manuel Kaucrs. Is the full counting function for gessel walks holonomic? Technical report, 

INRIA-Rocqucncourt, 2008. (in preparation). 
[3] Mireille Bousquct-Melou and Marko Pctkovsck. Linear recurrences with constant coefficients: the multivariate 

case. Discrete Mathematics, 225:51-75, 2000. 
[4] Frederic Chyzak and Bruno Salvy. Non-commutative elimination in Ore algebras proves multivariate identities. 

Journal of Symbolic Computation, 26:187-227, 1998. 
[5] Ira Gessel. personal communication, 2001. 

[6] Manuel Kauers and Doron Zeilberger. The quasi-holonomic ansatz and restricted lattice walks. Journal of 

Difference Equations and Applications, 2007. to appear. 
[7] Germain Kreweras. Sur une class de problemes lies au treillis des partitions d'entiers. Cahiers du B.U.R.O., 

6:5-105, 1965. 

[8] Marni Mishra. Classifying lattice walks restricted to the quarter plane. In Proceedings of FPSAC'07, 2007. 
[9] Neil J. A. Sloane. The On-Line Encyclopedia of Integer Sequences, http://research.att.com/~njas/sequences/. 
[10] Nobuki Takayama. An algorithm of constructing the integral of a module. In Proceedings of ISSAC'90, pages 
206-211, 1990. 

[11] Nobuki Takayama. Grobncr basis, integration and transcendental functions. In Proceedings of ISSAC'90, pages 
152-156, 1990. 

[12] Doron Zeilberger. A holonomic systems approach to special function identities. Journal of Computational and 
Applied Mathematics, 32:321-368, 1990. 

Manuel Kauers, Research Institute for Symbolic Computation, J. Kepler University Linz, Austria 
E-mail address: mkauersOrisc .uni-linz . ac . at 

Christoph Koutschan, Research Institute for Symbolic Computation, J. Kepler University Linz, Aus- 
tria. 

E-mail address: ckoutsch@risc .uni-linz . ac . at 

Doron Zeilberger, Mathematics Department, Rutgers University (New Brunswick), Piscataway, NJ, 
USA. 

E-mail address: zeilbergamath.rutgers.edu 



